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Abstract. In this paper we consider the Stratonovich reflected stochastic 
differential equation dXt = cr(Xt) o dWt + b(Xt)dt + dLt in a bounded domain 
O which satisfies conditions, introduced by Lions and Sznitman, which are 
specified below. Letting W^ 1 be the TV-dyadic piecewise linear interpolation 
of Wt what we show is that one can solve the reflected ordinary differential 
equation X? = a(X^)Wf + b(X^) + if and that the distribution of the 
pair (X^, Lf) converges weakly to that of (Xt, L t ). Hence, what we prove is 
a distributional version for reflected diffusions of the famous result of Wong 
and Zakai. 

Perhaps the most valuable contribution made by our procedure derives from 
the representation of Xj* 1 in terms of a projection of wt 1 . In particular, we ap- 
ply our result in hand to derive some geometric properties of coupled reflected 
Brownian motion in certain domains, especially those properties which have 
been used in recent work on the "hot spots" conjecture for special domain. 



1.1. Motivation. As is well known, ltd stochastic differential equations can be 
very misleading from a geometric standpoint. The classic example of this observa- 
tion is the Ito stochastic differential equation (SDE) 



where Wt is a 1-dimensional Brownian motion. If one makes the mistake of thinking 
that Ito differentials of Brownian motion behave like classical differentials, then 
one would predict the X(t) should live on the unit circle. On the other hand, Ito's 
formula, which is a quantitative statement of the extent to which they do not behave 
like classical differentials, says that d|A"(<)| 2 = \X(t)\ 2 dt, and so |X(i)| 2 = e*. 

To avoid the sort of misinterpretation to which Ito SDE's lead, it is convenient 
to replace Ito SDE's by their Stratonvich counterparts. When one does so, then the 
Wong-Zakai theorem [Mj shows that the solution to the SDE can be approximated 
by solutions to the ordinary differential equation (ODE) which one obtains by 
piecewise linearizing the Brownian paths. In this way, one can transfer to solutions 
of the SDE geometric properties which one knows for the solutions to the ODE's. 
The purpose of this paper is to carry out the analogous program for SDE's for 
diffusions which are reflected at the boundary of some region. This is not the 
first time that such a program has been attempted. For example, R. Petterson 
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proved in [7] a result of this sort under the assumption that the domain is convex. 
Unfortunately, convexity is too rigid a requirement for applications of the sort which 
appear in papers like [5] by Banuelos and Burdzy, and so it is important to replace 
convexity by a more general condition, like the one given in [5] by A. Sznitman and 
P.L. Lions. Finally, it should be mentioned that the article [5] by A. Kohatsu-Higa 
contains a very general, highly abstract approximation procedure which may be 
applicable to the situation here. 



1.2. Background for Reflected SDE's. We begin by recalling the (determinis- 
tic) Skorohod problem. 

Let O C R d be a domain and to each x £ dO assign a nonempty collection 
v(x) C S d_1 , to be thought of as the set of directions in which a path can be 
"pushed" when it hits x. Given a continuous path w : [0, oo) —> R d with w £ O, 
known as the "input," we say that a solution to the Skorohod problem for (0, f(x)) 
is a pair (x.,£.) consisting of a continuous path t £ [0, oo) i — > x t £ O and a 
continuous function of locally bounded variation t £ [0, oo) i — > £ t £ M. d such that 

(1) x t = w t +£ u \£\ t = f l ao (x s )d\£\ s , and 4= f v(x t )d\i\„ 

Jo Jo 

where \£\t denotes the total variation of £ t on the interval [0, i], and the third line 
is a shorthand way of saying that 

d£ t , , ,, „, 

6 v{x t ), d\£\t - a.e.. 



d\£\ 

When a unique solution exists for each input, we will call the map w. ~» (x.,£.) the 
Skorohod map and will denote it by F. Also, the path x. will be referred to as the 
"output." 

Throughout this paper we will take v{x) to be the collection of inward pointing 
proximal normal vectors 

(2) v[x) = {v£ S^ 1 : 3C> W £ 6 (x' - x) ■ v + C\x - x'\ 2 > 0}. 
Elementary algebra shows that 



(3) (a/-x)-v + C\x - x'\ 2 <0<=> x' - (x - — ? ) 



< 



2C 



which shows that, geometrically, u(x) is the collection of unit vectors based at x £ 
dO such that there exists an open ball touching the base of v but not intersecting 
O. 

The class of domains which we will consider was described by Lions and Sznitman 
in [5] . Namely, we will say that O is admissible if 

Definition 1.1. (1) Vx £ dO, v{x) ^ <j>, and there exists a Co > such that 
(x — x) ■ v + C Q \x - x'\ 2 > for all x' £ O, x £ dO, and v £ v(x). 
(2) There exists a function £ C 2 (]R d ;]R) and a > such that 
V(j>(x) ■ v > a > for all x £ dO and v £ v(x). 
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(3) There exist n > 1, A > 0, R > 0, a\, . . . , a n € § x , and x\, . . . , x n € <9<3 
such that 

n 

dO C (J B{ Xll R) and 

■i=i 

ieSOn B(x l ,2R) v ■ en > A > for all i/ G z/(x). 

In view of ([3]), Part 1 of Definition [Til can be seen as a sort of uniform exterior 
ball condition. More precisely, it says that not only can every point x € dO can be 
touched by an exterior ball but also that the exterior ball touching x can be scaled 
to have a uniformly large radius. In the convex analysis literature, the closure of a 
set O satisfying Part 1 of Definition [Li] is said to be uniformly prox-regular (See [5], 
especially Theorem 4.1, for more on the properties of uniformly prox-regular sets). 

Parts 2 and 3 of Definition |1.1| are regularity requirements on dO which ensure 
that the "normal vectors" don't fluctuate too wildly. In this connection, notice that 
Part 3 is implied by Part 2 when O is bounded. 

In their paper [5j, Lions and Sznitman show that for each w. € C([0,oo);R d ) 
there exists an almost surely unique solution (x.,£.) to the deterministic Skorohod 
problem when the domain O is admissible. The map T which takes w. to x. is 
called the deterministic Skorohod map. 

We turn next to the formulation of reflected diffusions in terms of a Skorohod 
problem for an SDE. Until further notice, we will be looking at Ito SDE's and will 
only reformulate them as Stratonovich SDE's when it is important to do so. 

Let O C R d an admissible domain, and let a : 6 — > Hom(M r ; R d ) and b : 6 — > 
WL d be uniformly Lipschitz continuous maps. Given an r-dimensional Brownian 
motion W. and x € O, a solution to (X , L.) to the reflected SDE Q is a continuous 
process {(Xj,L{) : t > 0} which is progressively measurable with respect to W. 
and satisfies the conditions that (X t , L t ) £ 6 x M. d and \L\ t < oo for alii > 0, and, 
almost surely, 



(4) 



X t =x + f a(X s )dW s + f b(X s )ds + L u 
Jo Jo 

\L\t= ! l 9 o{X 3 )d\L\ s , andL 4 = f u(X s )d\L[ 
Jo Jo 



where \L\ t denotes the total variation of L t by time i, and the third line is shorthand 

for arzft e K^t), d\L\ t -&.e.. 

Existence and uniqueness of solution to reflected SDE's was proved by H. Tanaka 
in [12] when O is convex. The extension of his result to admissible domains was 
made by Lions and Sznitman in 6J and Saisho in [TO] . We refer the reader to those 
papers for an overview of the subject. 



2. Equations with Reflection 

2.1. Properties of Solutions to Reflected ODE's. Suppose that O is abounded, 
admissible domain and that a : O — > Hom(R r ;R d ) is uniformly Lipschitz contin- 
uous. In this section we will show that, for each x € 6 and w. e C([0, oo);K d ) 
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there is precisely one solution (x. ,£.) to the reflected ODE 



x t =x + / a(x s )dw s + £ t , 
\t\t= / lao(x s )d\t\ s , and l t = / u(x s )d\£\ s 



where x. G C([0, oo); O) and £± : [0, oo) — > R d is a continuous function having finite 
variation \£\ t on [0, t] for all t > 0. In addition, we will give a geometrically appealing 
alternate description of this solution. Previously, existence and uniqueness results 
for variants of ^ are well known in the convex analysis literature. For example, 
see [3] for a recent such result as well as a good overview of other known results. 

Although the proofs of existence and uniqueness are implicit in the contents of 
other articles, we, mimicking the proof of Theorem 3.1 in [B], will prove them here. 
For this purpose, consider the map F w : C([0, oo);(5) — > C([0,oo);(5) given by 
Pw{y-) — r(^o + Jo o~(y s )dw s ), where T is the Skorohod map. We will show that F 
has a unique fixed point, and the key to doing so is contained in the next lemma. 

Lemma 2.1. For each T > there exists a C = C W (T) < oo such that for any 
pair of paths y. and y[ , 

(6) \F w {y.) t -F{y') t \ 2 < f e c{t ^\y T -y' T \ 2 dT for all t G [0, T] . 



Proof. Set z. — F(y.) and z[ — F(y'). Given T > 0, we will show that there is a 
C < oo such that 

\z t -z' t \ 2 <C^j\z T -z' T \ 2 dT + j\y T - y' T \ 2 dr^j , t € [0,T]. 

Once this is proved, the required estimate follows immediately from Gromwall's 
inequality. 

Let (f> be the function associated with O (see part 2 of Definition |1.1| . For any 
constant 7, we have that 

e -7[0(z t ) + 0(4)]rf( e 7[0(zt)+0«)]| Zt _ z £| 2 ) 

=2fa - z' t ) ■ [(o-(yt)dw t + d£ t ) - {a(y' t )dw t + <)] 

+ \z t - z' t \ 2 -f[V^(z t ) ■ (o-(yt)dwt + d£ t ) + W<t>{z' t ) ■ (a(y' t )dw t + dl' t )) 

= [(2(zt - z't) + y\zt - z' t \ 2 V^{z t )) ■ y{z t )} d\£\ t 
+ [(2(4 - z t ) + j\z t - ztfV^z't)) ■ v(z' t )]d\£'\t 

+ [2(z t - z' t ) ■ (a(y t ) - a(y' t )) + j\z t - z' t \ 2 (Vcp(z t )a(yt) + V<j>(z' t )a(y' t ))] dw t 

Taking 7 = — , we have that (cf. Part 1 of Definition 1.1 1 the first two terms 
are less than or equal to 0. Since a and V</> are Lipschitz continuous and ^ is 
bounded on finite intervals, we know that there exists aC = C W (T) < 00 such that 



<c( f \z T ~z' T \ 2 dr+ f \z T ~ z' T \\y T ~ x' T \dr + f \y r ~y' r \ 2 dr 
\Jo Jo Jo 



for t e [0, R]. Thus, because \z T — z' T \\y T — y' T \ < \\z T — z' T \ 2 + \\y T — y' T \ 2 , we get 
our estimate after replacing C by 2C. □ 



AN APPROX SCHEME FOR REFLECTED STOCHASTIC DIFFERENTIAL EQUATIONS 5 



Once we have Lemma 2. 1| one can apply a standard Picard iteration argument 
to show that F w has a unique fixed point and that this fixed point is the first 
component of the one and only pair (x t ,£t) which solves ([5]). 

We now want to describe a couple of important properties of the solution (x. ,£.). 

Lemma 2.2. Let (x.,£.) be the solution to |5|) for a given input w. and starting 
point xq G O. Then there exists a constant C , depending only on a, b, and O, such 
that 

d\x\ t < Cd\w\ t 

Proof. Set y t = xq + J q cr(x s )dw s . Then, x. = T(y.), and so it follows from Theorem 
2.2 in [B] that d\£\ t < d\y\ t . Since a is bounded on O, there exists a C < oo such 
that d\y\ t < Cd\w\ t , and therefore, because x t = yt + £t, we have that d\x\ t < 
d\y\t + d\£\t < C(d\w\t + d\w\ t ), from which the lemma follows immediately. □ 

We now introduce a more geometric representation of the equation |5]). For a 
closed set D C R d and z € R d , let = mf yeD \y — z\ denote the distance from 

z to D and denote by 

r g (z)^{,GR rf :liminf ^ (z + M =0} 

h\,0 h 

the tangent cone (a.k.a. the contingent cone) to D at z. Finally, let proj D (z) 
denote the (possibly multi-valued) projection of z onto D. 

The following is a version of a representation result which was introduced origi- 
nally in [1]. 

Theorem 2.3. LetO be a bounded, admissible set andw. a fixed, piecewise smooth 
input. If (x.,£.) is the unique solution to Q), then 

(7) x t =pToi T0(xt) (o-(x t )w t ), t-a.e. 

Conversely, given a solution x. to there exists an £. such that (x.,£.) is a 
solution to |5|). 

Remark 2.4. In general, the tangent cone Tu(z) is only closed and not necessarily 
convex. However, Part 1 of Definition |1.1| guarantees that Tq(z) is convex for all 
z E O (cf. Lemma 2.5 below) and so proj^(-) is single valued. 



In order to prove Theorem |2.3| we will need to introduce some concepts from 
convex analysis. For more information about these concepts and their properties, 
we refer the reader to the texts [9] and |13) . 

A non-empty set K C M. d is called a cone if v € K Xv e K for all A > 0. 
Given a cone K, we denote by K* its polar cone K* to be the set {w : v ■ w < 
0, Vw G K}. Next, for a given closed set D C M. d and a z £ D, we define the 
proximal normal cone to D at z to be the set 

N p D {z) = {v e R d : 3C> s.t. (y - z) ■ v < C\z - y\ 2 , Vy € D} 

and the Clarke tangent cone to D at z to be the set 

f D (z) = {v e R d : Vz„ € £> s.t. z n -> z, 3w„ e T D (z n ) s.t. u„ -> «}. 

Note that Trj(z) is always convex. 

We now present a lemma which records the properties of an admissible set O in 
terms of these concepts. 
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Lemma 2.5. Let O be admissible. Then 
(1 ) For each z £ dO , 



v £ N£(z) £ v{z) forv^O 



° \r 



(2) The graph of z — > Nq(z) is closed. That is, if Zi £ O, Vi £ N^Zi), Zi — > z 



and Vi — > v, then v £ Nq(z). 



(3) Tq(z) — Tq(z), and so it is convex for all z £ O. 

(4) N p {z) =T 6 {z)* for all z £ 6. 

Proof. 1. is immediate from our definitions. 

2. follows from 1. and Part 1 of Definition |1.1| Indeed, there exists a Co > 
such that for each i, 

(8) {z i -y)-v i + C (s \v i \\z i -y\ 2 > 0, Vy £ 6 

(note that when Zi € O, = and ([sj holds trivially). Taking i — > oo we see that 
(z — y) ■ v + Cq\v\\z — y\ 2 > for all y € (5, from which it follows that i; € N^(z). 

3. and 4. follow in a standard way from 2. See Chapter 4. of (T3] and Chapter 
6 of [5] (in particular Corollary 6.29) for the details. □ 



Using ideas from [3], we now prove Theorem 2.3 



Proof. (Proof of Theorem 2.3 1 First suppose (x.,£.) is a solution to From 
Theorem |2.2| and its proof, we see that x. and I. are locally Lipschitz and therefore 
that it = o~(xt)wt + it, t-a.e. Since Xt+h and Xt-h are in O, we have that 

(9) x t £ -T & {xt) n T (x t ), t-a.e., 

and, because TQ(x t ) is convex, i t is the projection of a(x t )vj t onto TQ(x t ) if and 
only if [ ajx Awt — it) ■ (v — i t ) < for all, t> € Tg(a;t). Note that by property 1. of 
Lemma 2.5 — ^ <E N^(x t ) (when xt € O this holds trivially), and so, by property 
4. of Lemma 2.5 and ([9]), we have that 

± t ■ it < o, ^ • i t > o =► i t • 4 = o. 

Therefore, using property 4. again, we have that 

(a(xt)w t ~ it) ■ (v - i t ) = -it ■ (v - x t ) = -It ■ v < 

as desired. 

Conversely, suppose x. is a solution to (JtJ) , and set £ t = J ' £ s — a(x s )w s ds. Then 
£o = and, since u is bounded, £. is a continuous function of locally bounded 
variation. Finally, because it is the projection of o~(xt)uit onto the convex set 
T^ixt), we have that 

-it ■ (v - it) = (o-(xt)wt - it) ■ (v - i t ) < 0, Vv £ Tb(x t ). 

Since it £ Tg^t) and Tg(a;t) is a convex cone, for each v £ Tg^t), a; t +w £ Tg(xt). 
Thus, by replacing v with v + x t in the inequality above, we find that —t t -v < for 
all v £ Tg(xt), and so — l t £ T^^xt)* = N^(x t ). Finally, by property 1. of Lemma 
2.5 this implies that (xt,it) is a solution to □ 
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3. Tightness of the Approximating Measures 

Let (C([0, oo); W), J 7 , W) be the standard r-dimensional Wiener space. That 
is, T is the Borel field for C([0,oo);R r ) and W is the standard Wiener measure. 
We will use W. to denote a generic Wiener path and Tt to denote the er-algebra 
generated by W. \ [Q,t]. Finally, for each positive integer N, let W N denote the 
TV-dyadic linear polygonalization of W.. That is, W^-n — Wm2 -JV and W N is 
linear on [m2~ N , (m + 1)2~ N ] for each m € N. 

Next, O C M. d will be a bounded, admissible domain, and b : O — > M. d and 
a : O — > Hom(K r ; K d )) will be uniformly Lipschitz continuous functions. Given 
a starting point xo £ O, for each W. , (X N ,L N ) will denote the solution to the 
reflected ODE Q with w t and a(x) replaced by, respectively, 

W *^j efxl and U^jj e Hom(R r xE:I d xl). 

{X£* : t > 0} and {L^ : t > 0} are then progressively measurable with respect 
to {F t : t > 0}, and we will use on the (X, L,W)-pathspace C([O,oo);0) x 
C([0, oo); R d ) x C([0, oo); R r ) to denote the distribution of the triple (Xf, if, Wf) 
under W. 

In first subsection, we show that the family {V N : N > 0} is tight on the 
(X,L, W / )-pathspace. In second subsection, we also develop some estimates which 
will needed for the next section. 

3.1. Tightness of the P^. By Kolmo gorov's Continuity Criterion, we will know 
that {V N : N > 0} is tight as soon as we prove that for each m G N and T > 
there exists a C m (T) < oo, which is independent of N, such that 

(10) E[\W t N -W s N \ 2m+1 ] < C m {T){t-sf m 

(11) E[LYf -Xff +1 ] < C m {T){t-sf m 

(12) E[\L?-L?r +1 ] < C m (T)(t-sf m . 



First note that (10 1 is an easy consequence of the equality E[|Wj — W s 

i + ln 



C m (t - s) 2 where C m = EjJ}Fi |2 



-N 



The proofs of (11) and (12) are a little more involved. 
Lemma 3.1. There is a C < oo such that for all s < i < s + 2" 
(13) \X t N -X?\<C\W t N -W s N \+C(t-s) 

Proof. When s and t lie in the same A^-dyadic interval, this follows more or less 
immediately from Theorem 1 2. 2 1 Namely, 

\X t N -X?\< \X N \ t - \X N \ S < C((\W N \ t - \W N \ S ) + (t - s)) 

= C(\W t N -W?\ + (t- s )) 1 

where the last equality comes from the fact that s and t lie in the same iV-dyadic 
interval. When they are in adjacent A^-dyadic intervals, one can reduce to the 
case when they are in the same A^-dyadic interval by an application of Minkowski's 
inequality. □ 
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It remains to handle s and t with t — s > 2~ N , and for this we will need the 
next two lemmas. Here, and elsewhere, \u\ is shorthand for the largest iV-dyadic 
number m2~ N dominated by u. That is, \u\ equals 2~ N times the integer part of 
2 N u. 

Lemma 3.2. For m > there exists a C m < oo such that for all s < t 

rt » 2 " 

rN xttN I jiwJVi 



(14) 
and 
(15) 



E 



E 



t 

(u-H) d|W w |„ 



<C m (i-s)^ 



< C m (i-s) 2 



Proof, li s < t lie in the same A^-dyadic interval we have that 

J* \W? - W{t } \d\W N \ u =4"|AWft | 2 jf'(u - \u\)du < 2 N \A W[ N sl \ 2 (t - s) 



(u-lu\)d\W"\ u =2*1^^1 / («-H)du<|A^j|(i- a ) 



and so 



(16) 



E 



E 



t s 2— 

(u- [u\)d\W N \ 



<C m (t- s) 



Cm . 2 



(*-«)< C m (t-s). 



Applying the Minkowski inequality, we see that the inequalities ( 16 1 continue to 
hold for general s < t. 



a 



i.i 



and set 7 



2C 



Lemma 3.3. Let <fi and a be as in Part 2 of Definition 
where Cq is the constant in Part 1 of that definition. Given s > ; there exist 
{J~t : t > 0} progressively measurable functions {Z TiS : t > s} and {V T . S : r > s} 
satisfying 

(17) \Z» S \ < C\X» - X?\, \Z^ S - Zlj < C\X» - X»\, and |t£| < C, 
with a constant C < 00, which is independent of s and N, such that 



(18) e^'W ~ *f| 2 < f Z».dW» + jf V u N s du for all t 



> s. 



Proof. Just as in the proof of Lemma |2.1 



< e i«x?) (2(Xf - X?) + 7 |Xf - X? | 2 V^(Xf )) a(X t N )dW t N 
+ e^(^) (2(X t N - X?) + 7 |Xf - Xf| 2 V0(Xf )) b(* t ")dt. 



from which ( 18 1 follows with 
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and 



VZ = (2(X% - X?) + \X» - xmsWX?)) ■ b(X»). 

Since V</>, 6, and a are Lipschitz continuous functions on the bounded domain O, 



it is clear how to choose the C in (17 1. 



□ 



We now prove ( 11 ) in the case that t — s > 2 by induction on m. Taking into 
account the fact that <j> is bounded, we can use ( 18 1 to derive the estimate 

2" 



E[|Xf - X?\ 2m+1 ] <C m E 



(19) 



CLE 



i \ 2" 



C m E 



t N 2" 



for some C m < oo. Because V^ s is bounded (see (17)), the third term is bounded 
by a constant times (t — s) 2 . For the first term we have that, for some constants 
C < oo, 



E 



JY 



< CE 



< CE 



< C(i-s) 2 



w»-w^\d\w»\ u 



CE 



(«- M)d|w"| u 



where the first inequality follows from (17), the second inequality from (13), and 
the third inequality from (14) and (15). Finally, for the second term we have that 

2" 



E 



< CE 



< C(t- 


s r 




< C(t- 


s r 




< C(t- 


s r 




< C(t- 


s r 





\Z^, s \ 2m du 



\x^-x»r d u 



(N 



[t-sf 



du 



du 



< C(t - s) 



where the first inequality is an application of Burkholder's inequality, the third 
inequality follows from (17), the fourth inequality is our induction hypothesis, and 
the fifth inequality follows from our assumption that t — s > 2~ N . Hence we will 
be done once we show that (11) holds when m = 0. But we can handle the base 
case by the same estimates as above, only now noting that the second term of ( 19 ) 
is in this case. 
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Finally, we must prove (12). Since 

we have that 

E[|Lf-Lff" 



dX t N = <r(Xf )dW t N + b{X?)dt + dLf 



] <CE[|Xf -Xf| 2 " ] + CE 



\{X»)dW^ 



+ CE 



\(X»)du 



We already know that the first term is bounded from above by C{t—s) 2 . Moreover, 
because b is bounded and < s < t < T, the third term is bounded above by a 
constant depending on T times (t — s) 2 . For the second term we have that 



E 



t N 2" 



o{XZ)dWt 



< CE 



(a(X^)-*(X^)dW» 



< CE 



< C E 



N 



X^\d\W"\ u 



CE 



+ C(t - 



\w. 



N 



W& } \d\W 



■E 



(u- \u\)d\W N \ u 



+ C{t-sf n < C(t-s) 2m 

where the second inequality follows is an application of Burkholder's inequality 
and the fact that a is bounded, the third inequality follows from (13), and the last 
inequality follows from ( 14 ) and ( [l5| . Putting these inequalities together we get 
pl. 



Given a ip : [0, oo) 



d , f3 € (0, 1], and t > s > 0, set 

\tp(u 2 ) - #Ui)J 
sup _ 

s<«i<«2<* («2-'"l)' 3 



As an immediate consequence of the estimates in (10), (11), and <[12j) combined 
with Kolmogorov's Continuity Criterion (cf. Theorem 3.1.4 in /citeStroockBook), 
we have the following theorem. 

Theorem 3.4. For each (5 < \, p e (1, oo), and T > 0, there exists a Kp^ p (T) < oo 
such that 

P{\\W N \\^ T] V \\X N \\ mTi V ||£ W ||^ [0 ,t] > R) < K , p (T)R-P for R>0. 

3.2. Controlling the Variation of L N . In general, the variation of a function 
cannot be controlled by its uniform norm. Thus, before we can apply the tightness 
result in the previous subsection to get the sort of result which we are seeking, 
we must give a separate argument which shows that the variation of L N can be 
estimated in terms of its uniform norm. To be precise, Theorem |3.5| says that the 
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variation of L N \ [0, t] can be estimated in terms of the uni form norm of L N \ [0, t] 
and the Holder norm of X N \ [0,t]. Hence, since Theorem 3.4 provides control on 
the Holder, and therefore the uniform, norms of the three processes W N , X N , and 
L N , our tightness result will sufficient for our purposes (cf. Theorems 4.1 below). 
In the following, and elsewhere, || tp\\ [t 1 ,t 2 ] = su Pre[ti,t 2 ] I^C 7 ")!- 

Theorem 3.5. For all < s <t, 

(20) \L N \ t - \L N \ S < C((t 8 )ir«||X"||i iM + 1)11^111.,*]. 

where R is the constant given in Part 3 of Definition 
Our proof follows the proof of Lemma 1.2 in [3]. 

Proof. Let 0\, . . . , O n den ote the open balls B{x\, 2i?), . . . , B(x n , 2R) appearing 
in Part 3 of Definition |l.l[ and choose an open set Oq so such that Oq Q O and 
O C O U ULi B(xi,R). Given x £ 6, let k(x) be the smallest 1 < k < n such 
that x £ B(xk,R), or otherwise let k(x) be 0. Next, set Co = s and define (, n for 
m > 1 inductively so that 

C m+ i = t A inf{r > Cm ■ X? i O k(X Nj}. 
Consider the time interval [CmiCm+l]- If Cm < t and k(X^ ) = 0, then L N \ 



L C 

lu+i - k " TI 

1, then (cf. Part 3 of Definition [□} 



[Cm, Cm+i] is constant and so \L n \q x — \L N \^ m =0. If Cm < t and k m = k(X^ ) > 



(if m+I " L?J ■ a km = f m+1 v(X?) ■ a km d\L N \ T > X(\L N \ Cm+1 - \L N \ Cm ). 

J Cm 

Hence, in either case, 

li^lu+x - l^lc™ < C|£f m+1 -£f m | < C||l"|| m . 

At the same time, if Cm+i < t and k(XgJ > 1, then \X^ - X^J > R and so 



R \X?^-X» 



< < \\X N \ 



(Cm+l — Cm) 4 (Cm+1 — Cn 

Thus if .M = sup{m : Cm+i < t}, then 



i[o,f]- 



■M v (t- S )||^||l 
^ < 1 + |{m : Cm+i < i and fc(*£) > l}| < 1 + ^ 4 4 ' [ ' 1 

which, in conjunction with the preceding, means that 
M-l 



\L% - \L% < £ (\L N \ Tm+1 - \L N \ T J + (\L% - \L"\ Tm ) 

m=0 

< (CM + 2)\\L N \\ M <C[(t- s)R- i \\X n \\i [Stt] +l]\\L N \\ [Stt] 



a 
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4. Associated Martingale 
and submartingale problems 

We know that the sequence of measures {P^ : N > 0} is on [X, L, 14 7 )-pathspace. 
Our eventual goal is to show that this sequence converges. Equivalently, we want 
to show that all limit points are the same. In this section we will show that every 
limit solves martingale and submartingale problems, and in the next section we will 
show that this fact is sufficient to check that convergence takes place. 

Up until now, we have needed only the assumptions that O is bounded and 
admissible, and a and b are Lipschitz continuous. However, starting now, we will 
be assuming that a € C 2 (0;Hom(IR r ;M d )). In addition, it will be convenient to 
make a change in our notation. Instead to writing the equation which determines 
(Xf,L?) (pathwise) as 

(21) dX? = a(X?)dW t N + b(X?)dt + dL? 1 X?=x , 
we will use the equivalent expression 

r 

(22) dX? = ]T V(X t N )d(W z N )t + V (X t N )dt + dLf , X» = x 

z=l 

where Vi is the ith column of the matrix a and Vq = b. At the same time, we 
introduce the vector fields V : 6 — > R d x R r given by V = ( * ) for 1 < i < r 



and Vq = I n ° ) , where {ei, . . . , e r } is the standard, orthonormal basis in R r . Then, 



x 




lAr -almost surely, 



(23) dYt=^Vi(X t )d(Wi)t + Vo(X t )dt, Y 



i=l 



where = I ^ J . In keeping with this notation, we use and Dy to 

denote the directional derivative operators on R d and M. d x M. r determined, respec- 
tively, by Vi and V,. Finally, for £ € K d , will denote the translation operator on 
C(R d x W; R) given by T^(x, y) = <p{x - y). 

Theorem 4.1. Let P be any limit point of the sequence {V N : N > 0}. Then for 
all h € C{;(R d x R r :R), 

(24) h(Y t ) - J* (^^[DlT Ls h](X s ,W s ) + [DyT Lt h](X s ,W s )\ ds 

is a F '-martingale relative to the filtration {B t '■ t > 0} generated by the paths in 
the (X, L, W)-pathspace. Also, for all f € C^(Pi d ; R) satisfying §£(x) > for every 
x G dO and v € v(x), 

(25) f(X t ) - f(x ) - J* (^J2d 2 v J(X s ) + D Vo f{X s )^ ds 



-suo-mar, 



tingale relative to the filtration {B t : t > 0} . 



We will begin with the proof of the martingale property for ( 24 ) , and, without loss 
in generality, we will do so under the assumption that h is smooth and compactly 
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supported. What we need to show is that for any limit point P, < s < t and 
bounded, continuous, immeasurable F : C([0,oo);IR d x R d x W) — !• [0,oo), 



(26) 



E 1 



h(Y t ) - h(Y s ) - / Ch{u)du F 







where we have used Ch(u) to denote the integrand in (24), and clearly it suffices to 
check this when s and t are M-dyadic rationals for some M <E N. Thus, it suffices 
to show that 



(27) 



h(Y t ) - h(Y s ) - / Ch{u)du F 







for M-dyadic s and t and bounded, immeasurable F £ C([0, co); R d xl d x R r ). 
For N > M, write 

2 N t-i 

h(Y t )-h(Y s )= fcQW)2-») - KY m2 - N ), 



and, for each term in the sum, use (23) to see that, P -almost surely. 



h(Y(m+l)2- N ) - h{ Y m2~ N ) 



Am+1)2 - > 

/ y j [DyT Lr h]{X T ,W r ){W i , m )dT 

Jm2- N r~f 



/ [DyT LT h](X T ,W T )dT, 
Jm2- N 



where W ittn = 2 N (W,((m + l)2~ N ) - Wi(m2- N )). 
Since 



2 N t r ( m +i)2- N ft 

V / [D v T L h](X Tl W T )dr= / [^ o T L >](X T ,M/ r )rfr, 

Jm2- N Js 



the second term on the right causes no problem. 
To handle the first term, note that 



[D v T L h] (X r , W T ) =[D v T Lm2 _ N h] (X r ,W r ) 

" E l [ D vTL„d Xk h] (X T , W T ) dh D 

Jm2- N 



k=l 



Since the second term on the right is dominated by a constant times |£|( m +i)2- 
\L\ m2 -N, we see that 



2^t / r (m+l)2-» \ 

E \J 2 n yPv^A ~ [DyT Lm2 _ N h]){X Tl W T ) dr\ W ilT . 



n (m+l)2 
Im2 

<C2-^E pN [\L\ t \\W\\ hm ] — >0 
as N — » oo. 



m=2 w s 
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Next, use (J22J) to see that 
[ D vTL m2 _ N h](X Tl W T ) = [DyT Lm2 _ N h](X m2 - N ,W m2 - N ) 

/ [D Vo DyT Lm2 _ N h] (X a ,W a ) da + Y. [ d ^ D v, T L m2 - N h] (X„ W a ) dL a 
+ E Wi,™ f T [ D v 3 D v T Lm2 _ N h] (X a ,W a ) da. 

Since the conditional P^-expected value of 

W iim [DyT Lm2 _ N h] (X m2 - N , W m2 - N ) 

given B s is zero, the first term on the right does not appear in the computation. 
Moreover, After integrating the second two terms over ^2-^,(771 + 1)2--^], multi- 
plying by Wi yTn , and summing from m = 2 N to m = 2 t, one can easily check that 
the absolute values of the resulting quantities have P^-expected values which tend 
to as N — > oo. 

Finally, again applying (22), one finds that 



/ T [ D vD v T Lm2 _ N h]{X a ,W a )da 

Jm2- N 



can be replaced by 

(r - m2- N ) [D v D v T Lm2 _ N h] (X m2 - N , W m2 -» ) 

plus terms which make no contributions in the limit as N — > oo. Hence, we are left 
with quantities of the form 

2 N t 

]T 2- 2N - 1 W J , m W l , m [DyDyT Lm2 _ N h](X m2 - N> W m2 - N ). 

m=2 N s 

Since the P^-conditional expected value of 2~~ 2N W j ^ m Wi :m is 2" <5~j,j, 



E 



£ r^Wj.nWi,™ [D Vj D Vi T Lm2 _ N h] (X, 



m2~ N > "m2 



W„ 



r)\F 



= — E 1 
2 



/ (m+l)2 JV 



2 ~ N E [ D vP v T Lm2 _ N h]{X m2 - Nl W m2 - N )\F 



m=2 N s 

which, as N —¥ oo, has that same limit as 

8., 



2 



j [Dy. Dy T Ls h) (X s , W s ) ds J F 



The proof of (25) is similar, but easier, and so we will skip the details. The only 
difference is that when we apply (22) to the difference /(^( m+ i)2-») ~ f(X m 2- N )> 
we throw away the dL T integral since, under our hypotheses, it is non-negative. 
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5. Convergence 

In this section we complete our program of proving the {P^ : N > 0} converges 
to the distribution of an appropriate Stratonovich reflected SDE. By the uniqueness 
result of Lions and Sznitman (Theorem 3.1 of [Hj) and the tightness which we proved 
in 3.1, the convergence will follow as soon as we show that every limit P is the 
distribution of that reflected SDE. 

Let P be any limit of {P^ : N > 0}. By Theorem we know that, for all 



h £ C'£{R d x 



(28) h(X t ~ L t ,W t ) - h(x , 0) - / Ch(s)ds is a P martingale, 

Jo 

relative to {B t : t > 0}, where 

£ M«) = \YX d %. t lM (X s ,W s ) + [DyT L h] (X S ,W S ). 



2 

i=l 



Using elementary stochastic calculus, it follows from (28) that {Wt ■ t > 0} is a 
P-Brownian motion relative to {Bt : t > 0} and that, P-almost surely, 

X t - x - ^ f i Y,[D Vi Vi\(X s ) + VoiX^ ds-L t = jT a(X s )dW s , 

which can be rewritten in Stratonovich form as 

(29) X t -X Q = Y / [ V z (X s )odW s + f V (X s )ds + L t . 

i=1 Jo Jo 

Thus, the only remaining question is whether {L t : t > 0} has the required prop- 
erties. That is, whether, P-almost surely, \L\ t < oo and L lo(X s ) d\L\ s = for all 
t > 0, and ^ e u(X t ) a.e. 

Since the local variation norm is a lower semi-continuous function of local uniform 
convergence, Theorem |3.5| tells us that, P-almost surely, L. has locally bounded 



variation. In fact, by combining that theorem with the estimates in Theorem 3.4 
one sees that, for all t>0, \L\ t has finite P-moments of all orders. 

In order to prove the other properties of L. we will use the second part of Theorem 
(O), which says that for every / G C2(R<i ; M) satisfying f£(x) > for all x e dO 
and v € v(x), 

(30) f{X t ) - [ Cf{X s )ds is a P sub-martingale 

Jo 

relative to {B t : t > 0}, where 

Cf(x)= l -j^Dlj{ X ) + D v J{x). 

i=l 



Now compare this to what one gets by applying Ito's formula to (29). Namely, his 
formula says that if f/ = / Q * Vf{X s ) ■ dL s then 

f(X t 



/ Cf(X s ) — £/ is a P-martingale. 
Jo 
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Thus, £. is P-almost surely non-decreasing. Starting from this observation and 
using the arguments in Lemmas 2.3 and 2.5 of one can prove the following 
lemma. 

Lemma 5.1. For f G C^(E. d ;M.), define as above. Then, P-almost surely, 
lo(X s ) d\£f\ s — 0. Moreover, if > for all x in an open set U and all 

v G v(x), then, P-almost surely, t L ljj(X s ) d£( is non- decreasing. 



Because • L. = it is obvious from the first part of Lemma 5.1 that 
J lo{X s ) d\L\ s = P-almost surely, and so all that we have to do is show that, 
P-almost surely, € v {Xt) a - e - To this end, let <p be the function in Part 2 of 
Definition |1.1| and define 

a f (x) = inf and b f (x) = sup |^ 

§£(*) v&v{x) d / v {x) 

for x G dO. 

Lemma 5.2. If {x n : n > 1} C dO, v n G v(x n ) for each n > 1, and (x n , v n ) — > 
(x,v) in dO x S^ -1 , then v G v{x). In particular, for each f G C^(M d ;]R), a? 
is lower semicontinuous and b^ is upper semicontinuous on dO. Furthermore, if 
(x, t) G dO x S^ -1 and there exists a (3 > such that V/(x) • I > fia' (x) for a set 
S of f G Cg(M. d ;M) with the property that {Vf{x) : f G S} is dense in R d , then 
iev[x). 

Proof. The initial assertion is an easy consequence of Parts 1 and 2 of Definition 
|Llj Next, suppose that x n — > x in dO. Because, by the first assertion, v(y) is 
compact for each y G dO, for each n > 1 there is a v n G v{x n ) such that o/(x„) = 
vi(^")-^" • Now choose a subsequence {x nm '■ m > 1} so that lim ^^,^ a? {x n ) = 
linim-^oo af(x nm ) and z/„ TO — > v in S^ -1 . Then G f(x) and so 

a f (x) < YJ} \ — < Uminf a f (x n ). 

\/<p\X) ■ V n->oo 

The same argument shows that is upper semicontinuous. 

Next, let (x,£) and f3 be as in the final assertion. Then, by Part 2 of Definition 
|l.l[ By taking / to be linear in a neighborhood of O, one sees that for every v G W 1 
there exists a v G v{x) such that v ■ I > • Hence, for each x' G O there is 

a v G f(a;) such that 

V (pyx) ■ v a 

which, by (|3|, means that I G v(x). □ 

Lemma 5.3. For each f G C?(R d ;R) ; P-almost surely d£: is absolutely continuous 
with respect to d£? and a? (X t ) < < b? [X t ) for d(f -almost every t > 0. 

Proof. First observe that / £: is linear. Now choose A > so that V(A</>— f)(x) ■ 
v > for all x G <90 and i/ G v{x). Then, £. A ^ - ^ = A£. — £r is P-almost surely 
non-decreasing, which proves that dt;? <C d^t and that < A P-almost surely. 



AN APPROX SCHEME FOR REFLECTED STOCHASTIC DIFFERENTIAL EQUATIONS 17 



The proof that, P-almost surely, a(t) = ^fy(t) lies between a? (X t ) and (X t ) 

for d£.^-almost every t > is a simple localization of the preceding. For example, to 
prove the lower bound, use the lower semicontinuity of a' to choose, for each n > 1, 
a finite cover of dO by open balls B(x/., n > ?"n)> 1 < k < k n such that x kn € <90, 
r„ < ^, and a f (y) > a f (x k , n ) - \ for all y G B(x k , n ,r n ) D <90. Then 

^-(y) > (a / (a; fc „) - ^)^-(y) for all 1 < k < k n , y G B{x kn ,r n ), and ^ G 
Now let /i be the Borel measure on C([0,oo);M d x R d x R r ) x [0, oo) determined 

by 

for all Borel subsets V of C([0, oo); R d x R d x M r ) and all a < b. Then, by Lemma 
5.l[ we can find a Borel measurable set A C C([0, oo); K d x K d x K r ) x [0, oo) whose 



complement has /i-measure and on which both 

X. G dO and l B(xfc B , rB) (X)(o'(a!* in ) - i) < l B(xh n>rn) (xA 

hold for all n > 1 and 1 < fc < fc„. Hence, again by the lower semicontinuity of a* , 
we see that > a(X). The proof of the upper bound is the same. □ 

Theorem 5.4. Let ¥ N be the distribution of (X N ,L N ) under Wiener measure. 
Then {P N : TV > 0} converges to the distribution P of the solution to the reflected 
stochastic differential equation (29). 

Proof. As we said earlier, everything comes down to showing that if P is a limit of 
{P N : N > 0} then, P-almost surely I = G v(X.) d\L\ .-almost everywhere. 
Thus, because, without loss in generality, we may assume that \l. \ = 1, the second 
part of Lemma 5.2 says that it suffices for us to show that, P-almost surely, there 
exist a /?. > such that Vf(X.) ■ I. > fi.a*(X.) d|L|.-a.e. for sufficiently many 
/'s. To this end, first note that, since is P-almost surely non-decreasing, (3. = 
V<f>(X.)-£. > d\L\.-a..e. P-almost surely. Second, because L. = J2f =1 £ Xi P-almost 
surely, we know that, P-almost surely, d\L\. <C d^ and that, for each / G C^(M. d ; M), 

(*) v/(x) ■i- = ^ = ^? v ^( x ) ' 1 * ^- a/ ( x ) d l L l- a - e - 

Finally, let D be a countable, dense subset of M d , and for each v G D choose 
/„ G C£(R d ;R) so that f v (x) = v ■ x in a neighborhood of 6. Then, P-almost 
surely, (*) holds simultaneously with f — f v for every v G D. □ 



Remark 5.5. In our derivation of Theorem (5.4) we used (30) to show that L. has 
the required properties. However, using the ideas in Lemma 1.3 of [6j, we could 
have based our proof on the fact that the approximating L N, s had these properties. 
Our choice of proof was dictated by two considerations. First, it seemed to us to be 
the simpler one. Second, and more important, it brings up an interesting question. 
Namely, does ( 30 1 by itself determine P? In [11] it was shown that ( 30 ) determines 
P when O has a smooth boundary and C is strictly elliptic, even if the coefficients 
are not smooth. Thus, the question is whether the same result holds when O is 
only admissible and the coefficients of C are smooth but may be degenerate. 
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6. Observations and Applications 

It should be noticed that although the approximating L N 's as well as limit L. 
have locally bounded variation, the we cannot replace our (X, L, W^)-pathspace 
with one in which the middle component is the space of continuous paths of locally 
bounded variation. The reason is that although L N will be absolutely continuous, 
L. will not. Indeed, consider reflected Brownian motion on the halfline [0, oo). 
In this case — sup 0<s<t [— ] is piecewise constant and therefore absolutely 
continuous. On the other hand, L t = sup 0<s<t [— W s ], which is the local time at 
of W. and as such is singular. 

The main application of our result that we consider is the following: Suppose that 
for each N, the paths satisfy a certain geometric property almost surely and 
the set S of paths which satisfy this geometric property is closed in C([0, oo); R d ). 
It then follows that the paths of X t also satisfy this geometric property almost 
surely since 



(31) 



P(5) > limsupP Ar (5) = 1 

N— >oo 



where, abusing notation, we use and P to denote the marginal distributions of 
P N and P on X-pathspace. That is, P^) = P N (Ax C([0, oo); R d ) x C([0, oo); R r )) 
and ¥{A) = V(A x C([0, oo); R d ) x C([0, oo); R r )). We conclude with several exam- 
ples of the sort of application which we have in mind. 



Example 6.1. In R 2 , let O be the rectangle [-1,1] x [0,2]. 
consider the Stratonovich reflected SDE 



Fix x a € O and 



where cr(x) 

(32) 
and 
(33) 



X2 
-Xi 



dX t =a(X t )odW t + dL t , X a = x 
Then 

If |x | > 1, \X t \ < \x \ for t > 



If |z | < 1, \X t \ = \x \ for t > 



a.s. 



-a.s. 



Proof. In view of (31), it suffices to prove that ( |32[) and (33) hold P N -a.s. The 
distribution of X. under P^ is, in view of Theorem |2.3[ the same as the distribution 
of X N under W, where X N solves the ODE 



X? = W o- ]To{xn {a{X?)W?) 



X, 



N 



x 



It is easy to check that Vx € O, a E R, x ■ proj T _ (a(x)a) is non-negative or non- 



positive according as \x\ > 1 or \x\ < 1. Hence, because, for each Wt, i (\X£ | > 

2X t N ■ proj T _ (<r(X t N )W t N ) dt-a.e., S and S for X N are obvious. Figure 

0°< x t ' I f I 

shows a sample path of X^ under W (to save space, we denote the "intended 

velocity" a(Xf)Wf by v t and the "actual velocity" proj T _ {xN) (a(Xf)W t N ) by 

a 



We next consider coupled reflected Brownian motion, for which we will need the 
following lemmas. 
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Figure 1. 



Lemma 6.2. Suppose O is bounded and admissible. Then O x O is bounded and 
admissible as well. Furthermore, for each (x,y) £ d(0 x O), the set of normal 
vectors v(x,y) defined by has the representation 

(34) v{x,y) = \[ x a ) : v x £ v{x),v y £ v{y),a\ + a\ = l,ai,a 2 > 



a 2 v. 



when (x, y) £ dO x dO, 



v{x, y) — | ^ g J : v x £ v{x) > , when (x, y) £ dO x O, 
and 

u(x, y) = | ( ® ~\ : v y £ v{y) \ , when (x, y) £ O x dO. 

Proof. The representation formulae are a straightforward consequence of the defi- 
nition of inward pointing unit proximal normal vectors in ^ . That O x O satisfies 
Part 1 of Definition |1.1| follows from the representation formulae and the fact that 
O satisfies Part 1 of Definition 11.11 

We next show that O x O satisfies Part 2 of Definition |1.1| Since O is bounded, 
<fi is bounded in O and so after adding a constant to (f> if necessary, we may assume 
that tj> > 1 in 6. 

Let y) = (j)(x)4>(y). Then for all (x, y) £ d(0 x £>), v £ u(x,y), we have, 
by our representation formulae, that 

V$(x, y) ■ v =ax<f>{y)V <f>{x) ■ v x + a 2 (j){x)V (f){y) ■ v y 
>ai4>(y)a + a 2 4>{x)a > a(a\ + a 2 ) > a 

(where (oi,a 2 ) = (1,0) and (0, 1) for the cases (x,y) £ (dO x O) U (O x dO)), and 
so Part 2 holds with the function &(x,y). Finally, as O x O is bounded, Part 3 
follows immediately from Part 2. □ 

Lemma 6.3. Let O be bounded and admissible. Then for (x, y) £ O x O, 

T &xo( x ,y) = T o( x ) x T o(v)- 
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Furthermore, 

mo\^ , A M - ( pr °jr O (»)(0 
P ro Jr 0x0 (»,„) v ^ j - v proj^fa) 

Proof. When D C M d is admissible, it follows from Part 3. of Lemma 2.5 that 



h\fl h 



(i.e. lim replaces liminf). Since O, and by Lemma |6.2| x O, are bounded and 
admissible, the hrst statement then follows immediately from the relation 



4x0 (( " )+M ,« ) ) = d%{x + hv) + d%{y + hw). 



x \ , / v 

y ) \w 

The second statement then follows from the first by a similar argument. □ 

6.1. Synchronously Coupled Reflected Brownian Motion. We now discuss 
synchronously coupled reflected Brownian motion. A d-dimensional synchronously 
coupled reflected Brownian motion is a 2<i-dimensional process Z t = (X t ,Y t ) in a 
product domain O x O which satisfies the reflected SDE 

dZ t = a(Z t )dWt + dL t , 

where 

' I 



a(z) 



I 



Note that, because a is constant, there is no difference between the Stratonovich 
and Ito versions of the above SDE. We will express this reflected SDE in a more 
convenient form as the pair of reflected SDEs 

dX t = dW t + dL u X = x and dY t = dW t + dM t , Y = y . 

We think of X t and Y t as being two d-dimensional processes which are driven by 
the same Brownian motion W t and which are constrained to lie in the same domain 
O. The two processes move in sync except for when one or the other is bumps 
against the boundary and gets nudged. 

We now consider the geometric properties of synchronously coupled reflected 
Brownian motion in two domains. Such properties were used to prove the "hot 
spots conjecture" for these domains (See [2] and [1] for more details). 

Example 6.4. Let O C K 2 be the obtuse triangle lying with its longest face on 
the horizontal axis, and denote its left and right acute angles by a and /?. Suppose 
x =/= y , and for x ^ y, let /.(x, y) = arg(y — x). Then, P-almost surely, 

(35) - < Z(x , yo) < a for all t either - /? < Z(X t ,Y t ) < a or X t = Y t . 

Proof. By pil, it suffices to show that S holds P^-a.s. Fix N and W t G fi. In 



view of Theorem 2.3 and Lemma 6.3 it will suffice to show that X t and Y t satisfy 



(351 where X? and Yf satisfy the ODE 



(36) . „ . 

Y t N =proj Ta(y/) (Wf), Y N = y 

It is straightforward to check that the functions X^ , Y t starting at Xj^ = 
xq, Yq = yo and defined inductively for t G {m2~ N , (to + 1)2 _JV ] by Xj = 
proj a (X^ 2 _ N + (t - m2- N )W t N ) and Y t N = proj (Y^ 2 _ N + (t - m2- N )W t N ) 



AN APPROX SCHEME FOR REFLECTED STOCHASTIC DIFFERENTIAL EQUATIONS 21 



satisfy (|36|. A simple geometric argument shows that if Z(x,y) <E [—(3, a] then 

^(P r °jo( a; )'P ro jo(y)) G hi 3 ."] or P TO io( x ) = P r °jo(2/)- Fr om this it follows by 
induction that and Y t N satisfy (35 1 as desired. Figure [2] shows a pair of sample 
paths X^ and Y t N in the interval m2~ N < t < (m + l)2 _Ar where we use v to 
denote the constant vector 2 Ar (M / ( m+1 ) 2 -« — W m2 -N). □ 




Figure 2. 



Example 6.5. (Proposition 2 in [T])We now consider synchronously coupled re- 
flected Brownian motion in a Lip domain. A lip domain is a domain in R 2 which 
is bounded below by a function fi(x) and above by another function fy^x) each 
of which is Lipschitz continuous with constant bounded by 1. The domains are so 
named because they look like a pair of lips (See Figure [3]) . 




Figure 3. 

Consider synchronously coupled reflected Brownian motion in a lip domain O 
where the defining functions fi(x) and f2{x) are smooth and have Lipschitz con- 
stants bounded by A < 1. Then O is a bounded admissible domain. Recall the 
definition of Z(x, y) from the previous example, and let x ,y £ K 2 be such that 
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%Q Vo and ^(xo,Vo) G [— f > f ]■ We have the following geometric property for the 
paths X t and Y t : 



(37) 



Vt, either Z(X t ,Y t ) e [-' 



or X = Y t 



-a.s. 



and Y t 



N 



Z(X t N 



satisfy (37 1 when X^ 1 and YJ 



Proof. In view of (|3T|) and Theorem |2.3| it suffices to show that for every € fl, 

N solve the ODE fl^. Let 8f ee 
It is enough 
tan -1 (A)] and 



or = YF. 



N ) or according to whether X^ ^ Y t N 
to show that, efa-almost everywhere, 8f < when 8^ G [f , f 
8^ > when —8^ € [f j § — tan" 1 (A)]. By symmetry it will suffice to prove the 
first statement. 

v t = WO)To(x t N )(vt), and v' t 



Let vt 



dt 



[8 



d 

— — tan 
dt 



N 



(Y t N 



X? 



{Y t N -X»)-R{v t 



l 

vt) 



P I0 iT o (Y t N )( v t)- We compute: 
_ (Yf-X t N ).R(v t -v' t ) 



\Y t N ~ 



*t N \ 2 



(Y»-Xn-R(v t -v>) 



\Y L 



N 



x?\ 2 



\Y t 



N 



where R = I ^ ^ J is the matrix which rotates vectors in M. 2 by 90° counter- 
clockwise. Suppose 6f G [f , f - tan" 1 (A)]. Then since the Lipschitz constants of 
fi and fi are strictly less than 1, Xf cannot be on the / 2 -boundary and Y t cannot 
be on the /i-boundary. For each t, it follows that either v t — v t or arg(i?(w t — v t )) G 
[71-— tan _1 (A), 7r+tan _1 (A)] and either v t — v' t or arg(R(vt — v' t )) G [n— tan _1 (A), n+ 
tan (A)]. And so each of the terms in the sum above is < 0. We depict in Figure 
[i] the case where X? G O and Y t N G 5C 




Figure 4. 

□ 

6.2. Mirror Coupled Reflected Brownian Motion. Our final example involves 
mirror coupled reflected Brownian motion. A <i-dimensional mirror coupled re- 
flected Brownian motion is a 2d-dimensional process Z t — (X t ,Y t ) in a product 
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domain O x O which satisfies the reflected SDE 



(38) 



dZ t = a(Z t )dW t + dL t 



where 



a(z) = a(x,y) 



1-2 



I 

\y-x\ 2 



defined up until the first time r that Z t hits the diagonal of O x O, at which point 
we stop our process (i.e. Z t = Z T for t > r). We will express this reflected SDE in 
a more convenient form as the pair of reflected SDEs 



(39) 



dX t =dW t + dL t , X Q 
dY t ={I - 2 



x 



(Y t - X t )(Y t - X t ) T 



)dW t + dM t , Y = y . 



We think of X t and Y t as being two d-dimensional processes which are "mirror cou- 
pled" with respect to the driving Brownian motion Wt and which are constrained 
to lie in the same domain O. That is, if you consider the hyperplane which perpen- 
dicularly bisects the line segment connecting X t and Y t to be a "mirror" , then the 
two processes move in such a way that they are mirror images of each other until 
either process bumps into the boundary and is nudged (which causes the mirror to 
shift). We refer the reader to the papers [5] and [T] for a more thorough overview. 

We will prove the same geometric property we considered for synchronously cou- 
pled reflected Brownian motion in Example |6.5[ but now for mirror coupled reflected 
Brownian motion. The point is that ( [38] ) can be viewed as a Stratonovich reflected 
SDE and so again it suffices to prove the geometric property for the approximating 
processes. 

We make this rigorous with the following lemma which shows that, off of the 



diagonal of O x O, the Stratonovich correction factor for (38) is 



Lemma 6.6. For t < r. 



(40) 




t (Y t -X t )(Y t -X t ) 
\Y-X f \ 2 



t)3 



In fact, 



(41) 



n (Y t -X t )(Y t -X t ) T ' 

i-i ) .(Wiw 



0, for each j 



Proof. It suffices to prove (41 ). Let Vi = (Y t — X t )i, where we have suppressed the 



dependence of Vi on t. An easy calculation shows that 



d(Vt,{W t ) 3 



z 2VjVi 



dt 
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and 



8 


( m 






d 


( 






d 


( v&i 







(Eft v k) 2 



Putting these together, we have that 



VjVj 

E k v h 




(E*^ 2 ) 2 
(E fc ^ 2 ) 2 

-IViVjVt 

(ET^fF 




From this, (41) immediately follows. 



□ 



We now prove a geometric property. 



Example 6.7. (Example 6.5 for mirror coupling) Let O be the same lip domain 



defined by smooth functions considered in Example 6.5 and consider the mirror 
coupled reflected Brownian motion starting from xq and yo where xq ^ j/o- Then 
( 37 1 holds where where X t and Y t are given by ([39]) . 



Proof. Let D e = {z = (x, y) G O x O : \x - y\ < e}_be the "e-diagonal" of O X O. 
Consider a sequence of smooth functions p k : O x O —> [0, 1] such that p(z) = on 
Di_ and p(z) = 1 off of Di . Let <T fc (z) = p k (z)a(z). Then cr fe G C 2 (<5 x (5). 

Let P fc be the measure on Z-pathspacc induced by the solutions to the reflected 
SDE 

dZ k = a k {Z k ) o dW t + dL k 

and define P k ' N to be the measures on Z-pathspace induced by solutions to the 
approximating reflected ODE 



(42) 



dZ^ N = a k {Z^ N )dWl 



N 



dL 



k,N 



Recall that the stopping time r corresponds to the first time X t equals Y t and define 
if{t : \X t - F t | < £}. Let S = {Z t G C([0, cx^); M 2d ) : -f < Z(X t ,Y t ) < 



Tk 



p2d\ 



]. :l ■ -\ and let S, \Z, ■ C([0, xi:„ , . , _^ >; . , / , . ,. 

Our goal is to show that P(5) = 1, where P is the measure induced on Z- 



■l<Z(X t ,Y t ) < f,Vt<r fe }. 



pathspace by (38). It is clear that the subsets S k decrease monotonically to S, and 
so it suffices to prove that P(S k ) — 1, Vfe. 
We first claim that V(S k ) = V k {S, 

a k on Di, it is clear that ¥(A) 



6.6 



k ) = ir~yo k ). mis if 
and the equality a 



and, in view of Lemma 

¥ k (A) for A G T Tk . So~we need only show that V k (S k ) = 1, and for this it will 
suffice to show that ¥ k ' N (S k ) = 1. We argue this as we did in Example |6.5| 
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Fix N and W t € fi and let 6 



k ' N = Z(X k ' N , Y k ' N ). By symmetry, it is enough 
to show that & k ' N < for Q k ' N G [f , | — tan _1 (A)] for almost every i < r^.. Let 



and 



w t 



I 



k.N 



X k < N )(Y k > N 



X. 



k,N\j ' 
t ) 



\Y 



k,N 



X 



fe.ATi 



Then, in view of Theorem 

proj To (Y t Jv ) (w t ) (recall that pk(X, 



Y t > =w t 



2.3 



and Lemma 



6.3 



k,N 



v t = proj To(X JV)(u t ) and 



'■",lf") = l for t<r k ). 



We compute: 



_ x N, k) . fl( - t _ - t) 



N,k 



lY N, k _ x N, k]2 



X^ K )-R{vt~v t ) , (y t 



w.fe 



+ 



\Y t 

N,k 



N.k 



X 



N.k, 



\Y 



N,k 



X 



N,k, 



X t N ' k ) ■ R(w t - w t ) 



\Y 



N.k 



X 



N,k\ 



4 yv t I 

The argument in the Proof [6~T] again shows that the first and third terms are non- 
positive. That the second term is non-positive follows from the fact that either 



w t or arg(w t - w t ) = ±6 



k.N 



a 
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